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A dynamically constrained coalescence model based on the phase space quantization and classical 
limit method was proposed to investigate the production of light nuclei (anti-nuclei) in non-single 
diffractive (NSD) pp collisions at ^/s—7 and 14 TeV. This calculation was based on the final hadronic 
state in the PYTHIA and PACIAE model simulations, the event sample consisted of 1.2x 10 8 events 
in both simulations. The PACIAE model calculated D yield of 6.247xl0~ 5 in NSD pp collisions 
at ^/s=7 TeV is well comparing with the ALICE rough datum of 5.456xl0 -5 . It indicated the 
reliability of proposed method in some extent. The yield, transverse momentum distribution, and 
rapidity distribution of the D, He, and ^H in NSD pp collisions at ^/s=7 and 14 TeV were predicted 
by PACIAE and PYTHIA model simulations. The yield resulted from PACIAE model simulations 
is larger than the one from PYTHIA model. This might reflect the role played by the parton and 
hadron rescatterings. 
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I. INTRODUCTION 

The nucleus-nucleus collisions at top RHIC energy 
produced an initial hot and dense matter (quark-gluon 
mater, QGM) and has been interpreted as a strongly cou- 
pled quark-gluon plasma (sQGP) [l|-0|. This is nearly a 
perfect liquid composed of quarks, antiquarks, and glu- 
ons bur is not a free gas-like quark-gluon plasma (fQGP) 
expected by theorists and experimentalists long time 
ago. 

We have proposed physical conjectures to theoreti- 
cally explore the property differences among the QGM 
formed in pp and/or nucleus-nucleus collisions at ultra- 
relativistic energies [5]. These conjectures are: (1) 
The fraction of gluon in yield and momentum increases 
with increasing reaction energy. (2) The quark flavors 
(u, d, s, ...) approach equilibrium with increasing reaction 
energy. (3) The antiparton to parton (u, d, s, ... quarks) 
ratio increases with increasing reaction energy. In short, 
the population of u,d,s,c,... quarks and the antiquark 
to quark ratio approach equilibrium (balance) and unity 
with increasing reaction energy, respectively. 

Recently STAR reported the observation of "an equi- 
librium in coordinate and momentum space populations 
of up, down, and strange quarks and antiquarks" by 
"The measured yields of \H (\H) and 3 He ( 3 He) are 
similar" or by the observation of %H /\H is close to 
3 He/ 3 He in Au+Au collisions at top RHIC energy @. 
Whether this STAR observation means the initial fire- 
ball created in this collision is really a fQGP may be 
debated. 

The investigation of anti-nuclei has great meaning in 
the nuclear and particle physics, the astrophysics, and 
even in the cosmology. Recently ALICE published their 
preliminary results of D production in pp collisions at 
y / s= 7 TeV [3] at nearly ten months later than STAR 
publishing their measurements for !^H production in 



Au+Au collisions at v /saT/v=200 GeV p. However, be- 
cause of quite low multiplicity, the study of anti-nuclei 
production is very hard both in experiment and theory. 

So far the report about the formation of QGM in early 
stage of pp collisions at RHIC energy is still absent in 
our knowledge. However, there were studies about the 
measurable of flow parameters and the elliptic flow sig- 
nature of the QGP phase transition in high multiplicity 
(energy) pp collisions |8H13|. So one could not rule out 
the QGM formation possibility in pp collisions at LHC 
or higher energies |14| . 

In this paper the PYTHIA model (H and PACIAE 
model [16J were used to calculate the final hadronic state 
in non-single diffractive (NSD) pp collisions at y / s=7 
and 14 TeV. That was followed by the calculations for 
the production of light anti-nuclei in the dynamically 
constrained coalescence model based on the phase space 
quantization and classical limit method. The PACIAE 
model calculated D yield of 6.247x 1CT 5 in NSD pp col- 
lisions at yfs=7 TeV is well comparing with the ALICE 
rough datum of 5.456xl0 -5 . This may indicate the re- 
liability of the proposed method in some extent. The 
parton and hadron rescattering effects were analyzed by 
comparing the PACIAE results with the PYTHIA one. 
It turned out that their role is un- negligible. 



II. MODELS 

PYTHIA is a model for high energy hadron-hadron 
(hh) collisions (151 . The parton and hadron cascade 
model, PACIAE[16}, is based on PYTHIA. In the 
PYTHIA model a hh collision is decomposed into parton- 
parton collisions. The hard parton-parton collision is de- 
scribed by lowest leading order perturbative QCD (LO- 
pQCD) parton-parton interactions with the modification 
of parton distribution function in a hadron. The soft 
parton-parton collision, non-perturbative phenomenon, 



is considered empirically. The initial- and final-state 
QCD radiations and multiparton interactions are consid- 
ered. So the consequence of a hh collision is a partonic 
multijet state composed of di-quarks (anti-diquarks) , 
quarks (antiquarks), and gluons, besides a few hadronic 
remnants. It is then followed by the string construction 
and fragmentation, one obtains a hadronic final state for 
a hh (pp) collision eventually. 

For the pp collisions the PACIAE model is different 
from the PYTHIA in follows: 

1. The string fragmentation is switched-off and the 
di-quarks (anti-diquarks) are broken randomly into 
quarks (antiquarks). So the consequence of a pp 
collision is a initial state of quarks, antiquarks, and 
gluons, besides a few hadronic remnants. This par- 
tonic initial state is regarded as the hot QCD mat- 
ter (QGM) formed in the relativistic pp collisions. 

2. The parton rescattering is introduced. In this stage 
the rescattering among partons in QGM is consid- 
ered by the 2 — ► 2 LO-pQCD parton-parton in- 
teraction cross sections [13] • However, a K fac- 
tor is introduced to consider the higher order and 
non-perturbative corrections. The effective strong 
coupling constant is assumed to be a s =0.47. A 
parton colour screen mass /^=0.63 GeV is intro- 
duced to avoid the divergence. Integrate the differ- 
ential cross sections above, the total cross section 
of the parton collision is obtained. Then the par- 
ton rescattering is simulated by the Monte Carlo 
method. 



formation of nuclei (hypernuclei) was not dynamically 
treated continuously but as a final state interaction sep- 
arately. The above analytical coalescence model relies 
upon the Wigner function [23] constructed according to 
the assumed wave function of light nuclei (anti-nuclei) 
[2l| . The statistical model relies upon the equilibrium 
and temperature assumptions as well as the restriction 
on projectile-like fragment in the intermediate energy 
heavy ion collisions. In this paper we proposed a dy- 
namically constrained coalescence model based on the 
phase space quantization and classical limit method to 
calculate the yield of light nuclei (anti-nuclei) after the 
transport model simulation directly. 

One knew in statistical mechanics [24| that, for a sys- 
tem with three dimensions in the state definition the set 
of microscopic states held in a volume element Ar cor- 
responds, in the limit of h — > 0, to a set of 



Ar 

7^ 



(1) 



quantum states. Of course, this correspondence is an 
approximation as long as h remains finite. In the Eq. [1] 
Ar reads 



Ar = AqAp 



(2) 



where q and p stand for the three position and three 
momentum, respectively. In other word, a h 3 volume 
element in the three dimensions phase space corresponds 
to a state of the system. 



3. The hadronization is proceeded after parton 
rescattering. The partonic matter can be 
hadronized by the Lund string fragmentation 
regime 1151 and/or phenomenological coalescence 
model (l6| . 

4. At last the hadron rescattering is added. In this 
stage the hadronic matter after hadronization suf- 
fers rescattering. It is dealt with by the usual 
two-body collision method [18| , until the hadronic 
freeze-out (the hh collision pair is exhausted). 

In short, the PACIAE model consists of the parton ini- 
tialization, parton evolution (rescattering), hadroniza- 
tion, and hadron evolution (rescattering) four stages. 

STAR Q has found in Au+Au collisions at top RHIC 
energy that "The measured \H/\H and 3 He/ 3 He ratios 
are consistent with the interpretation that the \H and 
\H are formed by coalescence of (A + p + n) and (A + 
p+n), respectively." Obviously, this simplest coalescence 
assumption is utilized only for ratio and not for yield. 

In the theoretical studies, the yield of light nuclei 
(anti-nuclei) was always calculated in two steps: In the 
first step the nucleons and hyperons were calculated by 
transport model. Then the light nuclei (anti-nuclei) 
were calculated by the analytical coalescence model [19l — 
|21| and/or the statistical model [22J. This meant the 



III. CALCULATIONS AND RESULTS 



Taking the system of ^-H as an example, a configura- 
tion consists of p, n, and A in a single event of the final 
hadronic state from transport model simulation could be 
expressed as 



CpnK(Qi,Q2,q3;Pi,P2,P3), 



(3) 



where subscripts 1, 2, and 3 were a shorthand for p, h, 
and A, respectively and q\ refers to the distance between 
p and the center-of-mass of p, n, and A, for instance. 
This configuration contributed an partial yield of 



1 if m inv < m ± Am, qi < Rq, 

<5i23 = { qi< -Ro, ?3 < Ro; 

otherwise: 



to the ^-H. In the above equation 



m mv = ](£i + E 2 + £ 3 ) 2 - (ft + p 2 + p 3 ) 



2ll/2 



(4) 



(5) 



is the invariant mass, mo and Rq stand for the rest mass 
and radius of the ^H, Ato refers to the allowed mass 

A ' 

uncertainty. The total yield of ^-H in a single event is 
the sum of above partial yield first over the Eq. [3] type 
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FIG. 1: (Color online) Light anti-nuclei transverse momentum distributions in the NSD pp collisions at y / i=7 and 14 TeV. 
The panels (a), (b), and (c) were for D (Am=0.0005 GeV), ^H (Am=0.005 GeV), and 3 He (Am=0.005 GeV), respectively, 
in the NSD pp collisions at y / s=7 TeV, the panels (d), (e), and (f) at y/s=14 TeV. The red solid and blue dashed histograms 
were calculated by final hadronic state in PACIAE and PYTHIA model simulations, respectively. 



TABLE I: Particle yield in NSD pp collisions at ^s=0.2 TeV. 





STAR 3 


PACIAE 


PYTHIA 


k+ 


0.140L0.010 


0.137 


0.125 


k~ 


0.137L0.009 


0.122 


0.115 


A 


0.0385L0.0036 


0.0382 


0.0309 


A 


0.0351L0.0033 


0.0381 


0.0312 



a The STAR data were taken from [26] 



configurations and then over the configuration types ob- 
tained by the combination among subscripts 1, 2, and 3. 
An average over events is required at last. 



The above summation was easily practised by follow- 
ing algorithm: One devised a three level loops (i, j, and 
k) over particle list in a single event of the final hadronic 
state from transport model simulation. If the particle 
code [H[ (KF code in PYTHIA model [Hj]) sum of the 
three particles (i, j, and k) was equal to -7446 (=-2212- 
2112-3122) and the condition of Eq. @]was satisfied then 
the ^-H yield increased by one unit. The loops was then 
proceeded by jumping to i loop, the regular loop pro- 
ceeding otherwise. 

In the PYTHIA and PACIAE simulations we assumed 
that the hyperons heavier than A will decay. The model 



parameters were fixed on the default values given in 
PYTHIA model, except the K factor and the parame- 
ters of parj(l), parj(2), and par] (3) (the later three were 
concerning in hyperon production [151 ]) were roughly fit- 
ted to the STAR data of k+k~,A, and A in NSD pp 
collisions at y/s= 0.2 TeV [26| as shown in Tab. Q] The 
fitted parameter values of 3 (default value is 1 or 1.5), 
0.15 (0.1), 0.38 (0.3), and 0.45 J0.4) were used to calcu- 
late the yield of D, 3 He, and ^-H in NSD pp collisions at 
y/s=7 and 14 TeV by the final haronic state in PACIAE 
and PYTHIA model simulations as shown in Tab. fll] 
One sees in the Tab. HT1 that: 

• The D yield of 6.247xl0~ 5 in NSD pp collisions 
at y/s—7 TeV in the PACIAE simulations seems 
reasonable comparing with the ALICE rough da- 
tum of ~ 5.456 x 10~ 5 (roughly estimated from 
Fig. 4 in [7] where integrating over p t and dividing 
by 350 M triggered events measured). The value 
in the PACIAE model larger than ALICE datum 
may attribute to the full rapidity phase space in 
the former (the same in all calculations in both 
PACIAE and PYTHIA models) bur |r?| < 0.9 in 
the later. 

• The PACIAE yield of anti-hadron increases with 
increasing reaction energy from 7 TeV to 14 TeV in 
a percentage of ~ 20%. It is less than the light anti- 
nuclei (~ 37% for D, - 50% for |i?, and - 70% for 
3 He). This might attribute to the available phase 
space increases with increasing reaction energy is 
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FIG. 2: (Color online) The same as Fig. [T] but for rapidity. 



TABLE II: Hadron and light nuclei (anti-nuclei) yields in 
NSD pp collisions at ^/s~7 and 14 TeV calculated by final 
hadronic state in the PACIAE and PYTHIA model simula- 
tions. 





PACIAE 


PYTHIA 




7 TeV 


14 TeV 


7 TeV 


14 TeV 


k+ 


4.563 


5.576 


3.802 


4.946 


k~ 


4.416 


5.331 


3.689 


4.778 


P 


4.152 


4.678 


3.491 


4.074 


V 


3.040 


3.588 


2.472 


3.078 


n 


4.094 


4.677 


3.397 


4.107 


n 


3.336 


3.938 


2.565 


3.335 


A 


1.648 


1.940 


1.285 


1.608 


A 


1.518 


1.769 


1.136 


1.386 


D" 


6.906E-05 


9.111E-05 


5.586E-05 


6.724E-05 


D° 


6.247E-05 


8.553E-05 4.852E-05 


6.048E-05 




5.456E-05 b 








\H C 


2.547E-07 


3.814E-07 


1.833E-07 9.677E-08 


W 


2.453E-07 


3.305E-07 


1.000E-07 


1.048E-07 


3 He c 


2.453E-07 


3.898E-07 


1.500E-07 


1.774E-07 


3 7U C 


2.642E-07 


4.407E-07 


1.417E-07 2.419E-07 



a Calculated with Am=0.0005 GeV. 

b Estimated from Fig. 4 in [7] . 

c Calculated with Am=0.005 GeV. 



stronger in the light anti-nuclei production rather 
than in the anti-hadron production. 

The yield resulted from PACIAE model simulation 
is larger than the one from PYTHIA model simula- 



tion at the same reaction energy. This reflects the 
role played by the parton and hadron rescatterings. 

Table IPLlI showed the feature of average transverse mo- 
mentum of the produced light nuclei (anti-nuclei) is not 
sensitive to the special piece of nuclei (anti-nuclei) and 
the reaction energy, like the case of hadron production. 
However, the (pt) seems to be ~ 1 Gev/c in light nu- 
clei (anti-nuclei) production instead of ~ 0.5 Gev/c in 
hadron production. This might mean the light nuclei 
(anti-nuclei) are generated more isotropically in the mo- 
mentum phase space. 

FigureCQgave the predicted D (Am=0.0005 GeV), |i? 

(Am=0.005 GeV), and 3 He (Am=0.005 GeV) trans- 
verse momentum distributions in NSD pp collisions at 
\fs= 7 and 14 TeV. The panal (a) was for D where 
the ALICE data pj were also given by the green circles. 
One sees in this panal that the dN/dpt shape in PA- 
CIAE and/or PYTHIA model simulations is similar to 
the one in general hadron production. The experimental 
Gaussian-like dN/dp t may be modified after "efficiency 
and annihilation corrections" [7J ■ The Fig. Q] show that 
the 1.2xl0 8 events are enough for dN/dp t distribution 
of D but not enough for ^-H and 3 He. This was also the 

° A 

reason using different Am value in the calculations for D 

and for ^H and 3 He. The larger fluctuation are shown 
A ° 

in panels (6),(c),(e), and (/). Globally speaking, the 
PACIAE dN/dpt distribution is not so much different 
form the PYTHIA one. 

We also predicted the rapidity distribution for the D, 
3 Se, and |iF in NSD pp collisions at ^/s=7 and 14 TeV 
in Fig. [5J In this figure one sees that the global fea- 
tures shown in dN/dy distribution are quite similar to 



the dN/dpt distribution in Fig. Q] so need not to repeat. 



TABLE III: Light nuclei (anti-nuclei) average transverse mo- 
mentum in NSD pp collisions at y/s=7 and 14 TeV calculated 
by final hadronic state in the PACIAE and PYTHIA model 
simulations. 





PACIAE 


PYTHIA 




7 TeV 


14 TeV 


7 TeV 


14 TeV 


D 


0.988 


0.993 


0.941 


0.952 


D 


0.999 
0.962° 


1.01 


0.955 
0.955 a 


0.981 


\H 


1.07 


1.05 


0.909 


1.01 


W 


1.14 


1.05 


1.09 


1.00 


3 He 


1.05 


1.07 


0.983 


1.04 


3 TH 


0.939 


1.10 


1.05 


1.05 



a Estimated according to the Fig. 4 in [7]. 



were employed to investigate the production of light 
nuclei (anti-nuclei) in NSD pp collisions at y/s=7 and 
14 TeV. We proposed a dynamically constrained coa- 
lescence model based on the phase space quantization 
and classical limit method to calculate the light nu- 
clei (anti-nuclei) by the final hadronic state in PACIAE 
(PYTHIA) model simulations. The calculated D yield 
of 6.247xlCT 5 in NSD pp collisions at ^1=1 TeV by 
PACIAE model final hadronic state is quite close to the 
ALICE rough datum of 5.456 xlCT 5 "7\. It turned out 
the reliability of proposed method in some extent. The 
light nuclei (anti-nuclei) yield, rapidity distribution, and 
transverse momentum distribution in NSD pp collisions 
at *fs=7 and 14 TeV were also predicted by PACIAE and 
PYTHIA model simulations. The yield resulted from the 
PACIAE model simulations is larger than the one from 
PYTHIA model which might reflect the role played by 
the parton and hadron rescatterings. 



IV. CONCLUSION 
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